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Oligomers of the organic semiconductor PTCDA are studied by means of helium nanodroplet 
isolation (HENDI) spectroscopy. In contrast to the monomer absorption spectrum, which exhibits 
clearly separated, very sharp absorption lines, it is found that the oligomer spectrum consists of three 
main peaks having an apparent width orders of magnitude larger than the width of the monomer 
lines. Using a simple theoretical model for the oligomer, in which a Frenkel exciton couples to 
internal vibrational modes of the monomers, these experimental findings are nicely reproduced. 
The three peaks present in the oligomer spectrum can already be obtained taking only one effective 
vibrational mode of the PTCDA molecule into account. The inclusion of more vibrational modes 
leads to quasi continuous spectra, resembling the broad oligomer spectra. 



I. INTRODUCTION 

Over the years there has been a lot of interest in the 
formation of finite molecular aggregates [ll-!6l|. In the in- 
vestigation of such aggregates optical spectroscopy plays 
an important role, since the aggregation manifests itself 
in (often drastic) changes of the absorption spectra 0- 
[Tsj . Usually in these studies the considered aggregates 
are embedded in a solvent at room temperature caus- 
ing considerable additional broadening of the individual 
vibronic absorption lines [l3 - [l6j . Although such opti- 
cal investigations already provide useful information e.g. 
on dimerisation constants, many details are not available 
due to the high temperatures and strong interactions of 
the aggregates with the environment (solvent). 

Using helium nanodroplet isolation spectroscopy it 
has become possible to study aggregation in an en- 
vironment (liquid helium) which interacts only weakly 
with the chromophores. As shown in Refs. [17i4l9| this 
technique allows to obtain monomer spectra of organic 
molecules with individually resolved vibronic lines. In 
this work we investigate absorption spectra of oligomers 
of the organic semiconductor PTCDA in helium nan- 
odroplets. PTCDA (3,4,9,10-perylene-tetracarboxylic- 
dianhydride, C24H8O6) is one of the most studied pery- 
lene derivatives. Its semiconducting properties and its 
ability to form hi ghly structured crystalline films on dif- 
ferent substrates [2^ makes this planar molecule inter- 
esting for applications ^2l] , such as organic light emit- 
ting diodes [22;, : thin-film transistors |24| - [27| or dye- 
sensitized solar cells (28l - [3]| . Most of the measurements 
on PTCDA concern thin films or crystals: highly ordered 
films of PTCDA have been grown and studied on poly- 
or sin gle crystalline substrates like mica and gold [33 |. 
silver |33| . graphene [3^ . aluminium, titanium, indium 
and tin |35| . Raman spectroscopy of epitaxial films (36| , 
electroabsorption [s^l , photoluminescence and Ultra- 
violet Photoelectron Spectroscopy (UPS) measurements 
on gold, silver and copper substrates [33| were also re- 
ported. 

Isolated PTCDA molecules have already been studied 



in different liquid solvents [40|, in solid Si02 matrix [4lj, 
on optical nanofibers [13] (sub-monolayer deposition) and 
with UPS in the gas phase [d^H^I. In all these cases there 
is a strong interaction of the PTCDA with an environ- 
ment. 

In the present study PTCDA molecules are embed- 
ded in helium nanodroplets where they are only very 
weakly perturbed, which makes them ideal systems to 
compare with theoretical model calculations. The iso- 
lation of organic molecules in helium nanodroplets has 
been established as a^j^rsatile technique for studying 
absorptionfsee Ref. 



45| and references therein) and 



emission [i^, properties at high resolution. Spec- 
troscopic measurements of single organic molecules and 
their derivatives with less than 1 cm^^ resolution in the 
vibronic spectra have already been reported, providing 
detailed information on vibrational structures for exam- 
ple for perylene [48,], PTCDA [l^, MePTCDI JjJ, Ph- 
thalocyanine and some polyacenes ^7|- Figure [T] 
shows a comparison of an absorption spectrum of single 
PTCDA molecules isolated in helium droplets (panel (b)) 
with the corresponding spectrum taken in an organic sol- 
vent (panel (a)), demonstrating the huge difference in the 
resolution of the vibronic transitions. 

Within the helium droplets several PTCDA molecules 
can aggregate to form oligomers. Such oligomers 
have been investigated by fiuorescence absorption spec- 
troscopy in order to study excitonic transitions 17l |. 
It was found that in contrast to the narrow discrete 
monomer vibronic absorption lines (cf. Fig. [IJb)) the 
oligomer spectrum, although in helium nanodroplets, has 
broad peaks, more like oligomer spectra in ordinary so- 
lution or on a surface. So far, there has been no assign- 
ment of the individual peaks and despite of the clean 
low temperature conditions, there has been no theoret- 
ical concept explaining the outcome and widths of the 
spectra. In Fig. [Ijc) a corresponding PTCDA oligomer 
spectrum is shown. While in the previous studies the 
oligomer spectra where measured only over a limited re- 
gion of excitation frequencies 17[, in the present work 
we could extend the covered range up to 24300 cm~^, 
revealing an additional peak at the high energy side. 



2 



20000 



21000 



22000 



23000 



24000 



PTCDA in an organic solvent 




(b) 



- PTCDA monomer 
in a helium droplet 
Gaussian convolution 




PTCDA oligomer 
in helium droplets 



20000 21000 22000 23000 24000 
Wavenumbers [cm"'] 



FIG. 1: (a) Measured absorption spectrum of PTCDA 
monomers in an organic solvent (DMSO) at room temper- 
ature taken from Ref. [131 • For better comparison the shown 
spectrum is shifted by 1800 cm~^ to higher wavenumbers. 
(b) Measured PTCDA monomer fluorescence excitation spec- 
trum in helium droplets [ItI] (solid line). The dashed curve 
is obtained by convoluting the measured spectrum with a 
Gaussian function having a width of 600 cm~^ to facilitate 
comparison with the spectrum shown in (a), (c) Measured 
PTCDA oligomer fluorescence excitation spectrum in helium 
nanodroplets 



In addition to extended experimental data, in this work 
we are able to reproduce the experimental findings using 
a simple Frenkel exciton model, that contains beside the 
dipole-dipole interaction between the monomers also the 
coupling of the electronic excitation to internal vibrations 
of the monomers. In particular, using this model, we will 
explain the appearance of the individual peaks of the 
oligomer spectrum and discuss the origin of their widths. 

We emphasize that the aim of this work is not a de- 
tailed fit of the measured oligomer spectrum but to gain 
basic insights into the underlying mechanisms. Therefore 
we adopt a minimalist model (for the PTCDA monomer 
as well as for the oligomer) , being aware that further re- 
finement could be necessary to obtain quantitative agree- 
ment between theory and experiment. 

The paper is organized as follows: In Section |TT] the ex- 
perimental setup is presented and details of the measured 
monomer and oligomer spectra are given. In Section IlIII 
we introduce the model of the monomer and the oligomer. 



Furthermore it is described how the absorption spectrum 
is calculated. In the following Section ITVl this theoretical 
model is applied to the calculation of the oligomer spec- 
tra. First the parameters needed to describe the PTCDA 
monomer are introduced. Since we cannot perform ex- 
act calculations of the oligomer spectra, including all the 
vibrational modes of PTCDA, in reasonable computing 
time, the concept of effective modes (EM) is introduced 
and illustrated. Using these EMs we then demonstrate 
the dependence of the oligomer spectra on the dipole- 
dipole interaction and also on the number of implemented 
EMs. With the insight gained from these considerations 
we then model the measured absorption spectrum. In 
Section |V] we summarize our findings and conclude with 
an outlook. 



II. EXPERIMENTAL SETUP AND FINDINGS 

The experimental setup is similar to the one presented 
in an earlier publication [l9| . Helium droplets are pro- 
duced by expanding highly pressurized helium gas (pu- 
rity 99.9999%, stagnation pressure 40 to 55 bars) into 
vacuum through a 10 fim nozzle kept at a temperature 
of 17 to 19 K by a closed-cycle refrigerator (Sumitomo 
RDK-408D). For these measurements we choose expan- 
sion parameters which lead to helium droplets with mean 
sizes of about 20 000 atoms In a second vacuum 

chamber, the helium droplets are doped with PTCDA 
molecules when flying through an oven placed about 
25 cm downstream from the nozzle aperture. For a fixed 
helium droplet size distribution, the number of PTCDA 
molecules picked-up by the droplets depends on the den- 
sity of PTCDA molecules in the oven given by its tem- 
perature and the known vapor pressure. The mixture of 
different oligomer sizes is given by the statistics of scat- 
tering events and follows a Poissonian distribution [50l |. 
At an oven temperature of 350 °C the droplets are prefer- 
entially doped with only single PTCDA molecules (max- 
imized monomer probability). For the oligomer spec- 
trum presented in Fig. H^c) the oven temperature was 
set to 390 °C. Poissonian statistics under these condi- 
tions reveal [l3| that approximately 8% of the droplets 
are undoped, 21% of the droplets are doped with a sin- 
gle PTCDA molecule, 26% picked up two molecules, 21% 
three, 13% four and 6% five molecules. In the given spec- 
trum (Fig. [ijc)) the very sharp monomer contributions 
have already been subtracted. Within the flight time 
to the detectors all embedded oligomers are cooled by 
evaporating helium atoms to the terminal droplet tem- 
perature of 380 mK [5l|. At these temperatures, only 
vibrational ground states are populated. Shifts of vibra- 
tional energies induced by the perturbation of the he- 
lium environment have been determined to be less than 
0.1% [13 . Shifts of electronic band origins of organic 
molecules embedded in helium droplets have been ob- 
served in the range of 10 - 100 wavenumbers j45| . 

About 30 cm further downstream, the doped droplet 
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beam is crossed by a pulsed laser beam for electronic ex- 
citation. The emitted Laser Induced Fluorescence (LIF) 
light is collected and focused on a photomultiplier placed 
perpendicularly to both the droplet beam and the laser 
beam. A second photomultiplier aside the imaged fluo- 
rescence spot monitors the stray light for efficient back- 
ground suppression. The droplet beam is chopped to 
allow for discrimination of PTCDA molecules coming di- 
rectly from the oven. In contrast to our former experi- 
ments we now used pulsed laser systems having the ca- 
pability to access a wider range of excitation frequencies. 
Two laser systems have been used. A dye-laser pumped 
by the 3'''^ harmonic (355 nm) of a Nd:YAG pulsed laser 
(Edgewave IS-IIIE) covers the green-blue range (500- 
370 nm). A second dye-laser system pumped by the 2"*^ 
harmonic (532 nm) of a pulsed Nd:YLF laser (Edgewave 
IS-IIE) provides output in the in the red-yellow region 
(590-500 nm). Typically pulses are used of ^ 10ns in 
widths at a repetition rate of 1 kHz having a spectral 
resolution of around 3 GHz. 



III. THE THEORETICAL MODEL 

In our model each PTCDA molecule is described by 
an electronic two level system (according to the 5*0 and 
state involved in the experiment). On each monomer 
the electronic excitation couples linearly to a set of har- 
monic vibrational modes. The oligomer is described by 
a Frenkel exciton Hamiltonian with transition dipole- 
dipole interaction between the monomers ^53 . j54i] . Other 
interactions (e.g. due to Van dcr Waals interaction, over- 
lap of electronic wavefunctions or interaction with the he- 
lium environment) are only considered via a global shift 
of the absorption spectra. 
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FIG. 2: Sketch of the harmonic potentials in ground and ex- 
cited electronic state of monomer n for one vibrational mode. 

n in the electronic ground state is then given by 

M 

and for the excited electronic state 
1 

= e» + 2 E (Pn, + <■ {Qn, - AQ„,f ) (2) 

where e„ denotes the energy difference between the min- 
ima of the upper and lower potential energy surface. 
We further introduce the dimensionless Huang-Rhys fac- 
tor Q 

EE ^(Ag„,)^ (3) 

which measures the strength of the coupling of electronic 
excitation of monomer n to vibrational mode j of this 
monomer. 



A. Hamiltonian of a single monomer 

First we consider a single monomer. Beside the elec- 
tronic ground state \4>^) we take one excited electronic 
state I (jf ) into account. In the Born-Oppenheimer ap- 
proximation we specialize to the particular case of har- 
monic vibrations of the same frequency in upper and 
lower potential energy surface. This allows for an efficient 
calculation of the oligomer spectrum (see appendix). 
Specifically, for each monomer we take M vibrational 
modes into account. The particular choice of the modes 
(or effective modes) will be discussed in Section ITVl For 
monomer n the coordinate belonging to the j-th mode 
is denoted by Qnj and P„j denotes the corresponding 
momentum operator. 

The vibrational frequency of mode j is denoted by w„j 
and the shift between the minima of the excited and 
ground state harmonic potential is denoted by AQ„j. 
The potential surfaces are sketched in Fig. [2] for the case 
of one mode. The vibrational Hamiltonian of monomer 



B. Hamiltonian of the aggregate 

One of our basic assumptions is that the N monomers 
forming the aggregate have fixed positions and orienta- 
tions w.r.t. each other. That means that we neglect vi- 
brational motion of the monomers w.r.t. each other and 
also conformational changes that could appear upon elec- 
tronic excitation. 

The electronic ground state of the aggregate is taken 
as the state 

|gei) = |0?)---|0^), (4) 

in which all monomers are in their electronic ground 
state. Since we consider linear absorption it is sufficient 
to take only states into account, in which at most one 
electronic excitation is on the aggregate, i.e. states 

JV 



4 



in which monomer n is electronicaUy excited and all other 
monomers are in their electronic ground state. 

To keep the model plain, we assume that in the aggre- 
gate the form of the monomeric potential energy surfaces 
is not changed. Furthermore we will not consider overlap 
of electronic wave functions of different monomers. As the 
important interaction responsible for the lineshape of the 
oligomer we take the transition dipole-dipole interaction 
between the monomers explicitly into account. In addi- 
tion we take energy shifts e.g. due to the chaiiged (Van- 
der-Waals) interaction between the monomers [5y,[53| via 
an overall shift of the calculated spectra into account. 
Since one expects that the individual energy shifts de- 
pend on the number N of monomers in the oligomer, we 
let also the included overall shift of the spectrum depend 
on N. 

The Hamiltonian in the electronic ground state of the 
aggregate is given by 



with 



N 



(6) 



N 



(10) 



These formulas are valid if all transition dipoles of the 
monomers are identical, which for simplicity we have used 
in the following calculations. 

Details of the numerical calculation can be found in the 
appendix. At this point we only want to note that ac- 
cording to Eq. (|9]) the spectrum is calculated by a Fourier 
transform of the correlation function c{t) = ( ^'o | ^'(0) 
with \ '^{t)) = e"''^"*/''! *o)- Practically this propaga- 
tion is carried out only to a certain finite time, which 
determines the resolution of the final spectrum. To guar- 
antee that c{t) decays to zero until the end of the propa- 
gation we multiply it with a Gaussian window function of 
width <Tt, which results in a convolution of the stick spec- 
trum with a Gaussian of the reciprocal width a — I/ctt- 
For the spectra shown in the following this width is indi- 
cated. 



We denote by Vnm the transition dipole-dipole interac- 
tion between monomer n and m, which is assumed to 
be independent of nuclear coordinates. The Hamiltonian 
in the one-exciton space defined by the states Eq. ([5]) is 
then given by 



N j N 

n— 1 \ m^n 

N 



(7) 



C. Calculation of the aggregate absorption 
spectrum 

Since the temperature of the molecules inside the he- 
lium nanodroplets is below 1 K, which is more than two 
orders of magnitude less than the temperature corre- 
sponding to the energy quantum of the vibrational mode 
with the lowest frequency (230 cm~^) of PTCDA, we as- 
sume that initially all vibrational modes of the PTCDA 
molecules are in their ground vibrational state. Thus, as 
the initial state of the aggregate we take the state 



<i>(t = 0)) = |5ci)|5vib), 



(8) 



where all monomers are in their electronic ground state 

and I .gvib ) = n«=i HjLi I "^"J = ^ ) '^i*^ I ""J = ^ ) de- 
noting the ground vibrational state of mode j of monomer 
n. We assume that the oligomers are isotropically ori- 
ented and we denote by the magnitude of the transi- 
tion dipoles of the monomers. The absorption strength 
for light with frequency v is then calculated from (ssl , Isoj 



A{iy) = N\nf Re 



dt 



-tH't/hi 



"fo) (9) 



IV. APPLICATION OF THE FORMALISM TO 
PTCDA OLIGOMERS IN HE NANODROPLETS 

A. The monomer parameters 

To apply the formalism of Section Hill we need to know 
the frequencies ujnj and the Huang-Rhys factors Xnj of 
the various modes of the PTGDA monomer. However, 
the numerical evaluation of Eq. ^ rapidly exceeds com- 
puter capabilities, even for the dimer, if many modes 
that couple to the electronic excitation are considered. 
For the calculations we therefore group several modes i 
with frequencies tJi that lie in a certain frequency interval 
together to form an effective mode (EM) j with a new 
Huang-Rhys factor 



(11) 



which is the sum of the Huang-Rhys factors Xi of the 
original modes i. This new EM j has a frequency 



.off 



1 



^3 ^ 



XiLOi 



(12) 



being the (Huang-Rhys factor weighted) average of the 
frequencies uji of the original modes i. This procedure has 
also been applied e.g. in Ref. [60]. However, one should 
keep in mind that introducing the EMs is an approxima- 
tion to the original spectral density. Of course, the more 
modes are combined to form an EM the more details of 
the original spectrum will be missing. We will discuss 
this point in more detail later in this section. 

Note also that the numerical effort for solving the 
oligomer problem not only depends on the number of 
EMs taken into account but also on the values of the 
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Wavenumber [1 /cm] Wavenumber [1 /cm] 

FIG. 3: Left column: spectral density of vibrational modes. 
Right column: corresponding calculated monomer spectra, a) 
All 31 modes of Ref. b) One EM combining the modes 

in the frequency interval 1000-1800 cm"^ of (a), c) 4 EMs 
combining all modes of (a), d) 6 EMs combining all modes of 
(a). The spectra in the right column are shown as sharp peaks 
and convoluted with a Gaussian of FWHM = 100 cm~^. 



Huang-Rhys factors. This is because the stronger the 
coupling of the electronic excitation to a vibrational 
mode is, the more vibrational states of this mode have 
to be taken into account in our basis for the numeri- 
cal calculation of the spectrum. Due to the fact that the 
Huang-Rhys factors of the EMs are the sum of individual 
Huang-Rhys (HR) factors of the original modes, the HR 
factors of the EMs will grow when reducing the number 
of EMs. For larger HR factors one needs a larger basis, 
which limits the number N of monomers of the oligomer 
that can be treated. 

When constructing the EMs we take the frequencies 
and Huang-Rhys factors extracted from the measured 
monomer spectrum Fig. [TJj, that are given in Ref. [l8l |. 
as starting point. Although it was noted in Ref. [3l that 
the peak intensities of the monomer spectrum do not 
exactly follow a Poissonian, as expected for harmonic 
modes, for the purpose of the present paper these de- 
viations are of no relevance. The frequencies ujj and 
corresponding Huang-Rhys factors Xj for all 31 modes 
found in Ref. are shown in Fig. [3^ and are listed in 
column A of Tab. HI Also shown is the corresponding 
monomer absorption spectrum obtained from the har- 
monic model (see Fig. [5^). For better visibility in the 
monomer spectrum the sharp peaks are also convoluted 
with a narrow Gaussian. 

An often adopted approach when considering 
oligomers in aqueous solution is to use only one EM 
describing the effect of the high energy vibrations 
\qi\ . corresponding to the progression of the broad 
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TABLE I: The four spectral densities used. Shown are the 
values of the mode frequencies iu (in cm~^) and the corre- 
sponding Huang-Rhys factors X. Column A: Values of all 31 
modes of Ref. [l8l | (shown also in Fig. |3^). Column B: Val- 
ues of the one EM (shown in Fig. ^jp). Column C: Values of 
the 4 EM (shown in Fig. [3j;). Column D: Values of the 6 EM 
(shown in Fig.[3Jl). The horizontal lines indicate which modes 
of column A are comprised in the respective EM in the other 
columns. 



maxima as seen in Fig. [Ija). In Fig. [31d we have illus- 
trated this approach, combining the high energy modes 
of the interval 1000-1800 cm^^ of Fig. to obtain one 
EM. Also shown is the resulting monomer spectrum (see 
Fig. [3f). Note the difference to the original spectrum 
Fig. [3^. We will use this model with only one EM to 
explain the three main peaks of the measured oligomer 
spectrum. 

To obtain a more realistic description of a monomer 
one would like to take as many EMs as possible. The 
number of EMs that can be taken into account in our 
numerical simulation strongly depends on the number 
of monomers within the oligomer. The spectrum of a 
dimer with 6-7 modes can be calculated within one hour 
on a standard PC. In a similar time, for a trimer we 
can only take 3 modes into account. In the calculation 
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of the spectra shown in the following we have used up 
to 6 modes for the dimer and 4 modes for the trimer 
(resulting in some increase in computer time). As noted 
above the grouping into effective modes is quite arbitrary. 
However, in the full spectral density Fig. [3^ one notices 
that several modes are densely packed together, like in 
the region around 1300 cm~^. Therefore it seems natural 
to put these modes together in one EM. Similar for the 
modes around 500 cm~^. The remaining high energy 
modes around 1600 cm~^ are also treated as one EM. 

Particular sets of four/six effective modes are shown in 
Fig. [3t/Fig. [3ji, which will be used in the following cal- 
culations. Of course the so constructed sets of EMs are 
only one of many possible choices, however quite reason- 
able ones. Note that the monomer absorption spectrum 
Fig. 13^ for the case of four effective modes resembles the 
full spectrum Fig. [5^ already quite well. 



B. Interpretation of the measured oligomer 
spectrum 

In this section we apply the foregoing model to explain 
the observed features of the measured oligomer spectrum 
Fig. [It. 

As described in Section |lT] the oligomer spectrum is 
expected to consist mainly of dimer and trimer contribu- 
tions. Since the computational effort strongly increases 
with the oligomer size and since we do not aim to describe 
small details of the measured spectra, we neglect in the 
following the contribution of oligomers consisting of more 
than three monomers. From the many possible geomet- 
rical arrangements we choose a particular simple one. 
We consider a stack of equidistantly spaced monomers 
with parallel transition dipole moments. Consequently 
we take the transition dipole-dipole interaction to be of 
the form Vnm — V/\n ~ m\'^ with a constant V. As 
pointed out in the previous sections all other interactions 
are only taken into account via an overall TV-dependent 
shift of the spectra. 



1. Model calculations 

To give an explanation of the appearance of the three 
peaks of the oligomer spectrum we first consider only one 
effective high energy mode per monomer with X and a; 
according to Fig. (i.e. X = 0.6 and oj — 1441 cm^^). 
Figure 2] shows dimer spectra (left column) and trimer 
spectra (right column) for various dipole-dipole interac- 
tion strength V (indicated on the individual figures). For 
small coupling V = ±200 cm^^ the spectra still look 
similar to the monomer spectrum (V — 0). The main 
difference is a small splitting of the high energy peaks, 
which lifts the degeneracy of the energy levels present 
for V — 0. After convoluting with a Gaussian one no- 
tices only a small change in the intensities of the resulting 
overall peaks upon going from V = QtoV = ±200 cm~^. 
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FIG. 4: Dimer (left column) and trimer (right column) spec- 
tra for different interaction strength V (in cm~^). Only the 
single vibrational mode of Fig. [3}3 is taken into account. The 
spectra are shown as sharp peaks and convoluted with a Gaus- 
sian of FWHM = 200 cm~^ (to facilitate comparison of the 
overall intensities). 



However, for increasing V this changes significantly. For 
V = 600 cm~^ we see in both cases, dimer and trimer, 
a peak structure (relative peak heights and distances be- 
tween the peaks) similar to that of the experimental spec- 
trum Fig. [TJ:. Thus with one effective high energy mode 
one can already understand the the three peaks of the 
experimental spectrum. 

If one takes now more effective modes into account one 
finds that a pronounced "broadening" of the peaks ap- 
pears. This is demonstrated in Fig. [S] where dimer and 
trimer spectra are shown for the same values of V as 
in Fig. 31 however, now including four effective modes 
per monomer. The used spectral density is that shown 
in Fig. [3J;; the values for Xj and ujj are given also in 
column C of Tab. [J Note that now in certain energy 
regions (e.g. around 23000 cm^^ for V — 600 cm~^) 
many peaks are closely spaced, forming a quasi contin- 
uum. This effect becomes especially pronounced in the 
high energy region of the spectra. The appearance of the 
many densely packed peaks leads to an effective broad- 
ening (but no significant change in position) of the three 
main peaks compared to the case when only one mode 
is taken into account. Thus to understand the broaden- 
ing of the aggregate spectrum it is of great importance 
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FIG. 5: Same as Fig. U but with four modes (according to 
Fig. [3j;) taken into account. 



to take as many internal modes as possible into account 
(this probably holds also true for spectra of PTCDA in 
other solvents, however there interaction with the solvent 
is expected to play an equally important role in broad- 
ening the spectra). 



2. Comparison to experiment 





21000 22000 

Wavenumber [1/cm] 



FIG. 6: a) Measured oligomer spectrum, b-d) Calculated 
oligomer spectra for V = 600 cm~^. b) Calculated dimer 
spectrum for the spectral density of Fig. [Sji, red-shifted by 
450 cm~^. c) Calculated trimer spectrum for the spectral 
density of Fig. [HJi, red-shifted by 550 cm~^. d) Sum of dimer 
spectrum (b) and trimer spectrum (c), where the two spectra 
are scaled according to the ratio 26/21 of dimer and trimer 
contribution to the experimental oligomer spectrum (a). The 
calculated spectra are shown as sharp peaks and in addition 
convoluted with a Gaussian of FWHM — 50 cm~^. 



A direct comparison between calculated spectra and 
the experimental one is given in Fig. [SI Here we have 
chosen V = 600 cm~^, for which the relative heights and 
distances of the three main peaks in the measured spec- 
trum (shown in Fig. [5^) are reproduced quite well. The 
dimer and trimer spectra (calculated including 6 EMs ac- 
cording to Fig. [5Ji for the dimer and 4 EMs according to 
Fig. [St for the trimer) are then shifted to the red as a 
whole so that the peak positions coincide with the mea- 
sured spectrum. We have used a slightly larger red-shift 
for the trimer than for the dimer (A = 550 cm"-'^ and 
A = 450 cm~^ respectively). As discussed in Section Hill 
such shifts are expected to be e.g. due to Van der Waals 
interactions. The so obtained dimer spectrum is shown 
Fig. [BJd, the trimer spectrum in Fig. Finally, in Fig- 
ure. [Gji a weighted sum of the dimer spectrum Fig.[6}D and 
trimer spectrum Fig. [SJ; is shown. The relative weights 
are 26/21 for dimer/trimer, taken to be consistent with 
the ratio of droplets containing two and three PTCDA 



molecules. To facilitate comparison with the measured 
oligomer spectrum we have convoluted the sharp peaks of 
the calculated oligomer spectrum with a narrow Gaussian 
of FWHM — 50 cm^^, much narrower than the apparent 
width of the three main peaks of the measured spectrum. 
Clearly, this artificial broadening has no physical reason. 
However, it facilitates the estimation of the overall inten- 
sity of the many sharp peaks. 

The foregoing considerations demonstrate that the ba- 
sic structure (relative heights, positions and widths of 
the apparent peaks) of the measured oligomer spectrum 
can be explained using a Frenkel exciton model, where 
intra- monomer vibrations are taken into account. 

In this work we have regarded the monomers to be in 
a fixed geometry. However, upon electronic excitation 
inter-monomer motion could be induced. The impor- 
tance of such motion has been demonstrated in Ref . 68:] . 
This will lead to an additional broadening of the Individ- 
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ual vibrational lines due to some low frequency motion 
(vibration/torsion/dissociation) of the monomers w.r.t. 
each other. 



V. CONCLUSIONS 

In this work we have presented a new extended LIF 
spectrum of PTCDA oligomers in helium nanodroplets. 
We have used a Frenkel exciton model including also 
intra-molecular vibrations to interpret the measured 
spectrum. Inclusion of one effective high-energy vibra- 
tional mode already allows to understand the appearance 
of the main features in the experimental data. Including 
more modes can explain at least part of the widths of the 
main peaks found in the experiment: While in the high 
energy region (above 21500 cm~^) of the theoretical spec- 
trum there is an already sufficiently large level density to 
explain the quasi-continuous measured spectrum, below 
21500 cm"-'^ the density of levels is still too low to match 
the experimental findings. Since we have restricted our 
calculations to dimers and trimers only, the inclusion of 
larger oligomers (which amount to roughly 1/4 of the 
measured oligomer spectrum) would give additional level 
density in this energy region. 

The used simple model apparently catches all the main 
characteristics of the spectrum despite of severe approxi- 
mations like e.g. fixed geometry and in particular not al- 
lowing charge transfer effects. In order to obtain a more 
detailed view of the type of oligomers present in the mea- 
sured spectra these additional effects should be taken into 
account. However, since the experimental data are not 
oligomer size selected, at the present stage it seems not 
to be needful to introduce a more sophisticated model 
risking over-interpreting the experimental data. Using 
hole-burning techniques it might be possible to obtain 
the individual contributions of dimers, trimers and higher 
aggregates and also to distinguish between different ge- 
ometrical arrangements. Then, sum rules [56j might be 
helpful to obtain information on the conformation of the 
different oligomers species. 

Nevertheless, the present study already reveals basic 
effects that will appear in the oligomer spectra. In par- 
ticular, it demonstrates that it is important to include as 
many internal vibrations as possible in the model calcu- 
lations. On the experimental side, vibrationally resolved 
spectra, as can be obtained by helium isolation spec- 
troscopy, appear to be instrumental as input for mod- 
eling and understanding the absorption of such organic 
complexes. 

Appendix A: Details of the numerical calculation 

1. Representation of the Hamiltonian 

We express the aggregate Hamiltonian using creation 
and annihilation operators for the vibrational modes. For 



the nuclear Hamiltonians of monomer n (see Eq. ([Ij and 
Eq. ©) one can write 

M 

and 

M 

(A2) 

where we have incorporated an energy X^jii f^nj^nj 
into e„. 

The electronic ground state Hamiltonian of the aggre- 
gate Eq. ([SJ is then given by 

(N M \ 
Y f^njol^anj I 5cl ) ( 5cl I (A3) 
«=lj=l J 

and for the Hamiltonian Eq. ([7]) in the one-exciton basis 
one has 

- Hel + i^int + i^vib (A4) 

with a purely electronic part 

N 

Hcl = Y i^riSnrn + Km) | 7r„ ) ( 7r„ | (A5) 

Ti,m— 1 

and a vibrational part given by 

N M 

Hvih ^YYj ^^njaljUnj. (A6) 
n=l j=l 

The coupling between the electronic excitation and the 
vibrations is contained in 

JV M 
Hint = - X I 7r„ ) ( 7r„ I X hWnj y/Xnj (a^^ + Qnj ) . ( A7) 
n=l j=l 

2. Calculation of the autocorrelation function 

To calculate the correlation function c{t) = 
(^-o |e-*^'*/'*| *o) , needed in Eq. © to obtain the 
absorption spectrum, we numerically solve the time- 
dependent Schrodinger equation 

thdt\^{t)}=H'\^{t)) (AS) 

with initial state | ^{t = 0) ) = | *o ) given by Eq. (fTO|) . 
To this end we expand the full state \'^{t)) in the basis 
given by the states 

ieM^K«)IW) (A9) 
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where the vibrational states | {a} ) are defined by 

N M 



w)-nn 



an 



(AlO) 



n=lj=l 



with the states \anj) being eigenstates of cinjCinj, i-e. 
'^lij'^njlanj) = anj\an]) with anj = 0,1,..., and the 
multi- index {a} contains all anj for all n and j. In 
this basis the matrix elements of the Hamiltonian H'^ 
are given by 

(N M 
fc=i j=i 

M c 

_ . V h'^}m . . . 



+ V„mS{a}{p}, 



where we have used the abbreviation 

N M 
n=l j=l 



(All) 



(A12) 



From Eq. (lAlip one sees that in the chosen basis the 
Hamiltonian becomes very sparse, because of the fac- 
tor S^a}{p} contained in each term. Finally, the time- 
dependent Schrodinger equation Eq. (jA8[) in this basis is 
given by 



For the numerical calculation we take only certain rel- 
evant states of the (infinite) basis, defined in Eq. (|A9p . 
into account. These states are chosen in the following 
way. First we restrict the maximum number each and- 
ean reach. In the state 1 7r„ ) , in which monomer n is 
electronically excited, wc take this number to be -^^j, 
so that anj < Z^^j. In a state | tt™ ) with n ^ m, in 
which monomer n is in its electronic ground state, we 
choose another maximum -, i.e. a„o < Zf ■. Further- 
more it has appeared to be advantageous to restrict the 
maximal number of vibrational states (and therewith the 
maximal vibrational energy) taken into account for a cer- 
tain j. We do this by setting a maximum Zj such that 

'^nj < Zj for the j-ih mode. The convergence of 
the calculated spectra is then tested by increasing one of 
the values -, Zf,, or Zj at a time. 

"■J "'J 



For example in the calculation of the trimer spectrum 
shown in Fig. |6t we took Zf = 2, Z| = 1, Zf = 2 
and Zf — 1 for the four used EMs given in column C of 
Table [T] (same values for all three monomers; the order of 
these values is the same as that of the respective mode 
frequencies in Table |l| and Z'^ = 11 for all modes and all 
monomers. Additionally in this calculation we made the 
restriction Zi = 3, Z2 = 2, Z3 = 3 and Z4 = 2 for the 
four modes. 

For the calculation of the autocorrelation function we 
use the method described in Ref. ^] to increase the ef- 
ficiency. 



et^ I ^(t) ; 



N M 

fc=i j=i 



M 



nj y Ckji 
. M 
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